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ABSTRACT 


The  term  M  out  of  N  sliding  window  detectors  refers  to  processes  that 
determine  whether  there  have  been  M  successes  in  a  sequence  of  N  binary  tri¬ 
als,  where  the  window  (of  N  trials  under  examination)  slides  along  a  possibly 
infinite  sequence  stopping  when  the  criterion  (of  M  successes  in  N  succes¬ 
sive  trials)  is  met.  They  are  frequently  used  to  model  the  operators  of  naval 
surveillance  systems  such  as  radar  and  sonar.  When  an  M  out  of  N  sliding 
window  detector  is  examining  a  sequence  of  trials  it  may  be  in  one  of  several 
states.  The  state  of  most  interest  is  the  accepting  state  into  which  it  enters 
(and  remains)  when  it  encounters  a  sequence  of  N  successive  trials  containing 
M  successes.  This  paper  describes  a  generalised  method  for  estimating  the 
probability  that  an  M  out  of  N  sliding  window  detector  is  in  its  accepting  state 
given  a  sequence  of  probabilities  representing  the  likelihood  of  success  on  each 
of  a  sequence  of  binary  trials. 
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Evaluating  the  state  probabilities  of  M  out  of  N  sliding 

window  detectors 


EXECUTIVE  SUMMARY 

The  term  M  out  of  N  sliding  window  detectors  refers  to  processes  that  determine 
whether  there  have  been  M  successes  in  a  sequence  of  N  binary  trials,  where  the  window 
(of  N  trials  under  examination)  slides  along  a  possibly  infinite  sequence  stopping  when  the 
criterion  (of  M  successes  in  N  successive  trials)  is  met.  They  are  frequently  used  to  model 
the  operators  of  naval  surveillance  systems  such  as  radar  and  sonar.  This  is  because  these 
systems  are  subject  to  false  alarms  and  the  operators  use  the  persistence  (scan-to-scan  in 
the  case  of  radar  and  ping-to-ping  in  the  case  of  sonar)  of  the  target  echoes  (as  one  cue)  in 
distinguishing  them  from  false  alarms.  As  the  echoes  are  intermittent  (especially  when  the 
target  is  at  extreme  ranges),  the  operators  are  quite  tolerant  of  missing  echoes  and  might, 
for  instance,  only  require  that  the  echo  appears  on  two  out  of  three  consecutive  scans/pings 
(in  which  case  they  would  be  performing  as  a  2  out  of  3  sliding  window  detector). 

When  an  M  out  of  N  sliding  window  detector  is  examining  a  sequence  of  trials  it 
may  be  in  one  of  several  states.  The  state  of  most  interest  is  the  accepting  state  which  it 
enters  when  it  encounters  a  sequence  of  N  successive  trials  containing  M  successes  and 
stays  there.  All  other  states  that  an  M  out  of  N  sliding  window  detector  may  enter 
while  examining  a  sequence  of  trials  are  known  as  non- accepting  states  and  this  includes 
the  initial  state  which  is  the  state  that  the  detector  is  in  before  any  trials  have  been 
evaluated. 

The  implementation  of  M  out  of  N  sliding  window  detectors  is  quite  simple  and  their 
incorporation  into  simulations  to  model  sensor  operators  is  straightforward.  Within  the 
simulation  the  appearances  (or  otherwise)  of  the  echoes  is  determined  and  the  result  is 
fed  to  a  simulation  of  the  required  sliding  window  detector  which  will  declare  a  detection 
when  the  criterion  is  met  and  enter  the  accepting  state. 

However,  in  studies  of  a  more  theoretical  nature,  what  is  known  about  the  series  of 
trials  (e.g.  radar  sweeps  or  sonar  pings)  is  not  the  result  (positive  or  negative)  of  each 
trial  but,  rather,  the  probability  that  the  trial  would  return  a  positive  result.  In  these 
cases,  what  is  required  is  an  estimate  of  the  probability  that  the  sliding  window  detector  is 
in  its  accepting  state.  A  generalised  (and  efficient)  method  for  solving  this  problem  (for 
arbitrary  values  of  M  and  N)  is  developed  in  this  paper  based  on  finite  state  automata. 
The  solution  is  simple  and  easily  implemented  in  modern  computer  languages. 
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1  Introduction 


The  term  M  out  of  N  sliding  window  detectors  refers  to  processes  that  determine 
whether  there  have  been  M  successes  in  a  sequence  of  N  binary  trials,  where  the  window 
(of  N  trials  under  examination)  slides  along  a  possibly  infinite  sequence  stopping  when 
the  criterion  (of  M  successes  in  N  successive  trials)  is  met.  They  are  frequently  used  to 
model  the  operators  of  naval  surveillance  systems  such  as  radar  and  sonar  (refs  1  and 
2).  This  is  because  these  systems  are  subject  to  false  alarms  and  the  operators  use  the 
persistence  (scan-to-scan  in  the  case  of  radar  and  ping-to-ping  in  the  case  of  sonar)  of 
the  target  echoes  (as  one  cue)  in  distinguishing  them  from  false  alarms.  As  the  echoes 
are  intermittent  (especially  when  the  target  is  at  extreme  ranges),  the  operators  are  quite 
tolerant  of  missing  echoes  and  might,  for  instance,  only  require  that  the  echo  appears  on 
two  out  of  three  consecutive  scans/pings  (in  which  case  they  would  be  performing  as  a  2 
out  of  3  sliding  window  detector). 

When  an  M  out  of  N  sliding  window  detector  is  examining  a  sequence  of  trials  it 
may  be  in  one  of  several  states.  The  state  of  most  interest  is  the  accepting  state  which  it 
enters  when  it  encounters  a  sequence  of  N  successive  trials  containing  M  successes  and 
stays  there.  All  other  states  that  an  M  out  of  N  sliding  window  detector  may  enter 
while  examining  a  sequence  of  trials  are  known  as  non-accepting  states  and  this  includes 
the  initial  state  which  is  the  state  that  the  detector  is  in  before  any  trials  have  been 
evaluated. 

As  will  be  shown  in  Section  2  the  implementation  of  M  out  of  N  sliding  window 
detectors  is  quite  simple  and  their  incorporation  into  simulations  to  model  sensor  operators 
is  straightforward.  Within  the  simulation  the  appearances  (or  otherwise)  of  the  echoes  is 
determined  and  the  result  is  fed  to  a  simulation  of  the  required  sliding  window  detector 
which  will  declare  a  detection  when  the  criterion  is  met  and  enter  the  accepting  state. 

However,  in  studies  of  a  more  theoretical  nature,  what  is  known  about  the  series  of 
trials  (e.g.  radar  sweeps  or  sonar  pings)  is  not  the  result  (positive  or  negative)  of  each 
trial  but,  rather,  the  probability  that  the  trial  would  return  a  positive  result.  In  these 
cases,  what  is  required  is  an  estimate  of  the  probability  that  the  sliding  window  detector  is 
in  its  accepting  state.  The  topic  of  this  paper  is  the  solution  of  this  problem.  The  key  to 
this  problem  of  evaluating  the  state  probabilities  of  M  out  of  N  sliding  window  detectors 
lies  in  examining  the  methods  for  implementing  them  and  to  this  end  the  next  section  will 
consider  their  implementation. 


2  Implementing  M  out  of  N  sliding  window 

detectors 


In  order  to  implement  an  M  out  of  N  sliding  window  detector  it  is  necessary  to 
maintain  a  history  of  the  last  N  —  1  trials.  If  the  result  of  each  trial  is  represented  by  a 
single  binary  digit,  with  1  representing  a  positive  trial  result  and  0  a  negative  one,  then 
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the  relevant  history  H  after  i  trials  is  given  by  the  string  (of  length  N  —  1)  of  binary  digits 

H  =  Ri-N+2Ri-N+3  *  •  -  Ri-iRi  (1) 

where  Rj  is  the  binary  digit  representing  the  result  of  the  jih  trial.  For  example,  if  N  =  4, 
and  the  result  of  the  last  three  ( N  —  1)  trials  was  positive,  negative  and  negative  (in  that 
order)  then  the  history  is  H  =  100.  In  order  to  determine  if  the  M  successful  trials  out 
of  N  consecutive  trials  criterion  has  been  satisfied  after  the  completion  of  a  trial  with  a 
positive  result  (given  that  the  criterion  has  not  already  been  satisfied)  it  is  necessary  to 
determine  the  number  of  digits  #  ( H )  with  the  value  1  in  the  history  H  and  if  it  equals 
M  —  1  the  criterion  is  satisfied.  If  the  criterion  was  not  satisfied  (after  a  successful  trial) 
it  is  necessary  to  update  the  history  H  in  preparation  for  the  next  trial  and  this  entails 
deleting  the  left  most  digit  from  H  and  adding  a  1  to  the  right.  For  example,  if  H  =  100 
before  a  successful  trial  that  fails  to  satisfy  the  criterion  then  it  becomes  H  —  001  after 
the  trial.  If  if  is  considered  to  be  the  binary  representation  of  an  integer  number  this  is 
equivalent  to  the  integer  arithmetic  operation  H  (H  x  2)  mod  2N~l  +  1  where  x  <=  y 
means  that  x  is  given  the  value  y  and  i  mod  j  is  the  modulus  of  i  with  respect  to  j. 

In  the  event  of  a  negative  result  to  a  trial  it  is  not  necessary  to  check  for  the  satisfaction 
of  the  M  out  of  N  criterion  as,  if  it  were  satisfied,  it  would  have  been  previously  determined. 
All  that  is  necessary  is  to  adjust  the  history  H  in  preparation  for  the  consideration  of  the 
next  trial.  This  necessitates  the  removal  of  the  left  most  digit  and  the  addition  of  a  0  to 
the  right  end.  If  if  is  considered  to  be  the  binary  representation  of  an  integer  number  this 
is  equivalent  to  the  integer  arithmetic  operation  H  <=  (H  x  2)  mod  2N~l. 


Algorithm  1  Conduct  binary  trials  until  M  successes  are  achieved  in  N  successive  trials 
accept  <=  false 
H  <=0 
repeat 

carry  out  binary  trial 

if  the  result  of  the  binary  trial  is  positive  then 

if  #  (H)  =  (M  —  1)  then  { M  out  of  N  criterion  is  satisfied} 
accept  4=  true 
else 

H  «=  (H  x  2)  mod  2N~l  +  1 
end  if 
else 

H  <=(Hx  2)  mod  2N~l 

end  if 

until  accept  =  true 


Algorithm  1  gathers  the  ideas  presented  above  to  provide  a  method  for  conducting  a 
series  of  binary  trials  until  the  criterion  of  M  successes  out  of  N  successive  trials  is  satisfied. 
Another  view  that  can  be  taken  of  this  process  is  that  it  is  a  finite  state  automaton  that 
transitions  from  one  state  to  another  dependent  on  the  result  of  a  binary  trial  until  it 
eventually  enters  the  accepting  state.  In  this  view  of  the  process,  each  possible  value  of 
the  history  H  can  be  considered  as  a  separate  state  of  the  automaton.  Let  Sm,n  be  the  set 
of  such  states  that  an  M  out  of  N  sliding  window  may  enter.  In  addition  to  the  accepting 
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state  there  will  be  one  state  in  this  set  for  each  integer  that  can  be  represented  by  N  —  1 
binary  digits  minus  those  integers  with  N  —  1  digits  that  contain  M  or  more  digits  with 
the  value  1  and  the  number  of  states  #Sm,n  in.  the  set  can  be  shown  (ref.  3)  to  be  given 
by 

N—M 

Tk  when  M  <  N  <  2M  -  1 

k— 0 


m,n  —  { 


M- 1 

Tk  when  N>2M-1 

k= 0 


(2) 


where 


Tk  =  ( 


)Af- 1 


)M- 1 


k-l  /w\l  / 

E(M\  i. 


N-M 
k 


when  k  =  0 


when  k  >  0 


(3) 


In  order  to  demonstrate  this  concept  consider  the  case  where  M  —  3  and  N  =  4.  There 
are  8  integers  (0,  1,  2,  3,  4,  5,  6  and  7)  that  can  be  represented  by  3  (i.e.  N  —  1)  binary 
digits  and  one  of  these  (7)  contains  3  (i.e.  M)  digits  with  the  value  1.  Let  the  label  of  each 
non-accepting  state  be  the  integer  interpretation  of  H  that  is  extant  when  the  system  is 
in  that  state  and  let  the  label  of  the  accepting  state  be  the  *  symbol.  Then,  using  state 
transitions  that  are  determined  by  the  processes  described  in  Algorithm  1,  the  finite  state 
automaton  that  implements  a  3  out  of  4  sliding  window  detector  is  described  by  the  state 
transition  table  shown  in  Table  1.  Alternatively,  it  is  described  by  the  state  transition 
diagram  shown  in  Figure  1  where  the  nodes  represent  the  various  states  and  the  edges 
represent  the  transitions.  It  is  this  view  of  the  system  as  a  finite  state  automaton  that 


From 

State 

To  State 

— ve 

-fve 

0 

0 

i 

1 

2 

3 

2 

4 

5 

3 

6 

* 

4 

0 

1 

5 

2 

★ 

6 

4 

★ 

★ 

★ 

★ 

Table  1:  State  transition  table  for  3  out  of  4  sliding  window  detector 


is  the  basis  for  the  evaluation  of  M  out  of  N  sliding  window  detectors  when  the  actual 
results  of  each  binary  trial  are  not  known  but,  instead,  the  probabilities  that  each  of  a 
sequence  of  binary  trials  is  successful  is  known. 
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0,1 


Figure  1:  State  transition  diagram  for  3  out  of  4  sliding  window  detector 

3  Evaluating  M  out  of  N  sliding  window 

detectors 

If  the  notation  Ps  ( fc )  is  used  to  represent  the  probability  that  the  process  is  in  state 
s  after  the  fcth  trial  then  the  probability  that  the  process  is  in  the  accepting  state  after 
the  fcth  trial  is  given  by  P*  (fc).  In  order  for  the  system  to  be  in  state  s  after  the  fcth  trial 
either  the  fcth  trial  must  have  had  a  positive  result  and  the  system  was  in  a  state  where 
it  transitions  to  the  state  s  after  a  positive  result  or  the  kth  trial  had  a  negative  result 
and  the  system  was  in  a  state  where  it  transitions  to  state  s  on  a  negative  result  and, 
therefore, 

Ps  (k)  =  Pkp?  (k  - 1)  +  (i  -  Pfc)  pr  (k  - 1)  (4) 

where  p*  is  the  probability  that  the  fcth  trial  has  a  positive  result,  P+  (A:  —  1)  is  the 
probability  that  the  system  is  in  a  state,  after  the  ( k  —  l)th  trial,  that  will  transition  to 
state  s  if  the  fcth  trial  has  a  positive  result  and  P~  (A;  —  1)  is  the  probability  that  the 
system  is  in  a  state,  after  the  ( k  —  l)th  trial,  that  will  transition  to  state  s  if  the  fcth  trial 
has  a  negative  result. 

For  example,  in  the  case  where  M  =  3  and  N  =  4  the  positive  transition  probabilities 
can  be  determined  to  be 


P0+(k~  1) 

=  0 

(5) 

P+(fc-l) 

=  P0(fc-1)  +  P4(fc-1) 

(6) 

P2+(fc-l) 

=  0 

(7) 

P3+(k~  1) 

=  Pi{k-1) 

(8) 

Pt(k~  1) 

=  0 

(9) 

P?{k~l) 

=  P2(fc-  1) 

(10) 

P+(fc-l) 

=  0 

(11) 

P*+(k~  1) 

=  p3(fc  — l)  +  p5(fc  — l)+p6(fc  — l)+p*(fc-l) 

(12) 

from  Table  1  and/or  Figure  1.  Similarly  the  negative  transition  probabilities  can  be 
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determined  to  be 


P0~(k-  1) 

=  P0(fc-1)  +  P4(fc-1 ) 

(13) 

Pf(fc-l) 

-  0 

(14) 

P2-(*-l) 

=  Pi(fc  — 1)+P5(fc-1) 

(15) 

P3~(k-1) 

=  0 

(16) 

Pi(k~  1) 

=  p2(fc-l)+p6(fc-l) 

(17) 

P5~(k~l) 

=  0 

(18) 

Pf(k-1) 

=  P3(k- 1) 

(19) 

P*(k~l) 

=  P*(k~  1) 

(20) 

and  this,  in  turn,  allows  the  full  set  of  state  probability  equations  for  a  3  out  of  4  sliding 
window  detector  to  be  determined 


Po(fc)  =  (l-p*)(Po(fc-l)+P4(fc-l))  (21) 

Pi(fc)  =  pfc(P0(fc-l)  +  P4(fc-l))  (22) 

P2(k)  =  (1  -pk)  {Pi  (£  -  1)  +  Ps  (k  -  1))  (23) 

P3(k)  =  pfcPi(fc-l)  (24) 

P4(fc)  =  (1  —  Pk)  {P2  (fe  —  1)  +  Pe  {k  —  1))  (25) 

P5(k)  =  pkP2  (k  -  1)  (26) 

P6(fc)  =  (l-Pjfe)P3(fc-l)  (27) 

P*(fc)  =  P*  (A;  —  1)  +  pk  (P3  (&  —  1)  +  P5  (A;  —  1)  +  P6  (A;  —  1))  (28) 


which  can  be  used  to  evaluate  the  detector  when  started  with  the  following  initial  values: 


Ps  (0)  = 


when  s  =  0 
otherwise. 


As  the  detector  remains  in  the  accepting  state  once  it  has  entered  that  state  the  probability 
p *  ( k )  that  it  enters  the  accepting  state  as  a  result  of  the  fcth  trial  is  given  by 


P*(k) 

P*(fc)-P*(fc-  1) 


when  0  <  k  <  1 
otherwise. 


(30) 


and  this  is  equivalent  to  the  probability  density  function  for  P*  (k).  It  can  then  be  used  to 
determine  the  most  probable  detection  range(s).  The  set  of  equations  needed  to  evaluate 
M  out  of  N  sliding  window  detectors  with  other  values  for  M  and  N  can  be  obtained  by 
similar  processes  to  that  just  described. 

Although  they  will  correctly  evaluate  M  out  of  N  sliding  window  detectors,  the  meth¬ 
ods  just  outlined  suffer  from  one  severe  disadvantage.  Every  time  it  is  required  to  evaluate 
a  sliding  window  detector  for  values  of  M  and  N  that  have  not  been  previously  encoun¬ 
tered,  a  new  set  of  equations  must  be  generated.  What  is  required  is  a  more  general 
method  that  can  automatically  accommodate  many  combinations  of  M  and  IV  and  this 
will  be  developed  in  the  next  section. 
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4  A  General  Method  of  Evaluation 


In  order  to  generalise  the  methodology  described  in  the  previous  section  each  possible 
state  of  a  system  will  be  represented  by  an  ordered  hexad  (if,  P,  H +,  if-,  P+,  P~),  where 
if  is  the  states  identifying  number,  P  is  the  probability  that  the  system  is  this  state,  if + 
is  the  identity  of  the  state  reached  from  this  state  after  a  positive  trial  result,  H~  is  the 
identity  of  the  state  reached  from  this  state  after  a  negative  trial  result  and  P+  and  P“ 
are  provided  to  hold  intermediate  values  during  the  evaluation  process.  The  set  S m,n  is 
the  set  of  such  hexads  necessary  to  represent  all  of  the  states  of  an  M  out  of  N  sliding 
window  detector  and  the  notation  Sm,n  [j]  is  used  to  refer  to  the  member  of  S m,n  that 
has  if  =  j.  Where  necessary,  dot  notation  will  be  used  to  signify  components  of  a  hexad 
(e.g.  the  notation  Sm>n  \j\  -P  refers  to  the  P  component  of  the  hexad  Sm,n  L?])* 


Algorithm  2  Evaluate  M  out  of  N  sliding  window  detector 
§m,a  [0]  .P  <=  1 

for  each  hexad  h  in  S m,n  where  h.H  ^  0  do 

h.P  <=  0 
end  for 
loop 

determine  the  probability  p  of  a  positive  result  on  the  current  trial 

for  each  hexad  h  in  S  m,tv  do 
h.P+  <=  0 
h.P'  4-  0 

end  for 

for  each  hexad  h  in  S  m_,n  do 
S m,n  [h.H+]  ,P+  4=  S m,n  [h.if+]  P+  +  h.P 
S m[n  [h.H~]  P-  4=  S m,n  [h.H~]  .P~  +  h.P 

end  for 

for  each  hexad  h  in  S  m,n  do 
h.P  <=  p  (fc.P+)  +  (1  -  p)  (h.P-) 

end  for 

wait  for  next  trial 

end  loop 


Algorithm  2  illustrates  a  process  by  which  the  set  of  hexads  Sm,n  uiay  be  used  to 
evaluate  an  M  out  of  N  sliding  window  detector.  After  initialising  each  hexad’s  state 
probability  P  to  the  appropriate  values,  it  enters  an  infinite  loop  and  evaluates  each  of 
the  state  probabilities  once  the  probability  of  a  positive  result  from  the  current  trial  is 
known.  At  any  time  that  the  system  is  waiting  for  the  next  trial  it  may  be  interrogated 
to  determine  the  various  state  probabilities.  In  particular,  the  value  of  Sa/?at  [*]  .P  will  be 
the  probability  that  the  system  has  successfully  met  the  M  out  of  N  criterion  based  on 
the  trials  processed  so  far  and  is,  therefore,  in  the  accepting  state. 


Algorithm  3  Create  the  set  of  hexads  S m,n 
S Af,N  =  {(*5  0,  ★,  ★,  0, 0)} 

Add  a  hexad  with  H  =  0  to  Sm,n  —  Algorithm  4 
Sjif  ,jv  [0]  .P  <=  1 
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In  order  for  Algorithm  2  to  be  applied  it  is  first  necessary  to  generate  the  set  of  hexads 
Sm,n  and  Algorithm  1  provides  the  key  to  solving  this  problem.  Based  on  the  information 
in  that  algorithm,  the  recursive  algorithm,  Algorithm  4,  and  Algorithm  3  can  be  used  to 
construct  the  set  Sm,n-  Algorithm  3  initialises  S m7n  to  contain  a  single  hexad  representing 
the  accepting  state  and  then  uses  Algorithm  4  to  add  a  hexad  for  the  initial  state  H  =  0 
to  Sm,n . 


Algorithm  4  Add  a  hexad  with  H  =  i  to  S m,n 
create  new  hexad  h 
h.H  <=  i 
h.P<=  0 

if  #  (i)  =  (M  —  1)  then 
h.H+  ^  * 

else 

h.H+  <=  (i  x  2)  mod  2JV_1  +  1 

end  if 

h.H~  <£=  (i  x  2)  mod  2N~X 

S  M  ,JV  <=  Sm,n  U  {h} 

if  a  hexad  with  H  =  h.H^  does  not  exist  in  S  m,n  then 
Add  a  hexad  with  H  —  h.H+  to  Sm,n  —  Algorithm  4 
end  if 

if  a  hexad  with  H  =  h.H~  does  not  exist  in  Sm,n  then 
Add  a  hexad  with  H  =  h.H~  to  Sm,n  —  Algorithm  4 
end  if 


The  recursive  nature  of  Algorithm  4  then  ensures  that  hexads  for  all  of  the  other  states 
that  an  M  out  of  N  detector  requires  are  added  to  S m,n  with  their  state  probabilities  set 
to  0.  Finally  Algorithm  3  sets  the  state  probability  for  state  0  to  1  ensuring  that  §m,n  is 
in  a  valid  configuration  (with  the  sum  of  the  state  probabilities  equal  to  1).  The  number 
m,n  of  hexads  in  the  set  S m,n  generated  by  these  algorithms  is  given  by 

#Sm7n  =  #S  m,n  (31) 

(see  Equation  2)  and  may  be  quite  large. 

5  A  General  and  Optimal  Method  of  Evaluation 

From  Algorithm  2,  it  can  be  readily  determined  that  the  time  taken  to  evaluate  the 
system  each  time  a  trial  is  processed  will  be  proportional  to  the  number  #Sm,n  of  hexads 
in  Sm,n  (which  is  given  by  Equation  2  and  can  be  very  large  depending  on  the  values  of 
M  and  N).  Therefore,  if  it  is  possible  to  reduce  the  number  of  hexads  in  Sa^at,  without 
altering  the  value  obtained  for  S m,n  [*]  -P  after  each  trial  is  processed,  then  the  evaluation 
process  can  be  made  more  efficient.  That  this  is  possible,  at  least  in  the  case  where  M  —  3 
and  N  =  4,  can  be  determined  by  examining  Table  1  and  Figure  1  and  observing  that 
two  states  are  effectively  equivalent  if  their  transition  states  are  identical.  In  the  case 
where  M  =  3  and  N  =  4  the  states  0  and  4  both  transition  to  state  0  on  a  negative  trial 
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result  and  to  state  1  on  a  positive  trial  result.  They  are,  therefore,  effectively  equivalent 
and  can  be  merged  into  a  single  state  to  form  a  new  finite  state  automaton,  described  in 
Table  2  and  Figure  2,  that  is  also  capable  of  determining  when  there  have  been  3  trials 


Prom 

State 

To  State 

— ve 

+ve 

0 

0 

i 

1 

2 

3 

2 

0 

5 

3 

6 

★ 

5 

2 

* 

6 

0 

* 

* 

★ 

* 

Table  2:  State  transition  table  for  an  optimal  3  out  of  4  sliding  window  detector 


Figure  2:  State  transition  diagram  for  an  optimal  3  out  of  4  sliding  window  detector 


with  positive  results  out  of  4  successive  trials.  Additionally,  an  equivalent  set  of  equations 
to  Equations  21  to  28 


P0(k) 

=  ( l-pk)(Po(k-l)  +  P2(k-l)  +  P6{k-l )) 

(32) 

Pi  (k) 

=  P*(Po(*~l)) 

(33) 

P2(k) 

=  ( l-pk)(Pi(k-l)  +  P5(k-l )) 

(34) 

Ps(k) 

=  PkPi(k~  1) 

(35) 

Ps(k) 

=  PkP2(k-  1) 

(36) 

Pe(k) 

=  (l-j>*)P3(fc-l) 

(37) 

P*(k) 

=  P*  (k  ~  1)  +  Pk  (P3  (k~l)+  Pa  (k  - 1)  +  P6  (k  -  1)) 

(38) 

could  be  generated  from  this  new  automaton  in  order  to  evaluate  3  out  of  4  sliding  window 
detectors  (if  desired). 

In  general,  the  process  of  reducing  a  finite  state  automaton  to  its  optimal  form  will 
be  more  complex  than  this  example  and  several  iterations  may  be  required.  Take,  for 
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example,  the  case  where  M  —  5  and  N  =  6.  The  finite  state  automaton  generated  using 
the  methods  described  in  Section  2  is  given  in  Table  3(a)  and  contains  31  states.  Among 
these  states  there  are  11  pairs  of  equivalent  states  (16  and  0,  17  and  1,  18  and  2,  19  and 
3,  20  and  4,  21  and  5,  22  and  6,  24  and  8,  25  and  9,  26  and  10,  and  28  and  12)  and  if 
each  of  these  pairs  is  merged  into  a  single  state  (taking  the  smaller  of  the  state  labels 
as  the  label  for  the  new  state)  then  the  automaton  described  in  Table  3(b)  is  obtained. 
Note  that,  even  though  all  the  equivalent  pairs  of  states  were  merged  into  single  states, 
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Table  3:  State  transition  tables  for  5  out  of  6  sliding  window  detectors 
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there  are  still  four  pairs  of  equivalent  states  (namely  8  and  0,  9  and  1,  10  and  2,  and  12 
and  4)  in  this  new  automaton.  If  the  same  process  is  repeated  on  this  new  automaton 
the  automaton  described  in  Table  3(c)  is  obtained  and  this  automaton  also  contains  an 
equivalent  pair  (4  and  0)  of  states.  After  one  more  iteration  of  the  process  the  automaton 
described  in  Table  3(d)  is  obtained  and  (finally)  this  is  an  optimal  automaton  with  no 
pairs  of  equivalent  states. 

If  an  iterative  process  based  on  the  above  discussion  is  applied  to  the  set  of  hexads  S m,n 
produced  by  Algorithm  3  then  the  optimal  set  S'M  ^  (that  which  represents  an  automata 
for  implementing  M  out  of  N  sliding  window  detector  in  the  minimal  number  of  states) 
can  be  generated.  SrMN  can  the  be  substituted  for  S m,n  hi  Algorithm  2  to  provide  an 
optimal  algorithm  for  evaluating  M  out  of  N  sliding  window  detectors.  The  number  of 
hexads  in  this  optimal  set  is  given  by 


which  in  most  instances  is  much  smaller  than  that  given  by  Equation  31  for  How¬ 

ever,  the  time  taken  to  generate  this  optimal  set  (using  this  process)  will  be  proportional 
to  #S m,n  rather  than  N  and  what  is  required  are  versions  of  Algorithms  3  and  4  that 

generate  the  optimal  set  directly. 


Algorithm  5  Create  the  optimal  set  of  hexads  N 

S'm,N  =  {(*7  O7  *7  *7  0)} 

Add  a  hexad  with  H  =  0  to  the  optimal  set  S'M7V  -  Algorithm  6 
S'm.jv  [0]  P  <=  1 


Algorithms  5  is  a  modified  version  of  Algorithm  3  which  sidesteps  the  problem  by 
shifting  the  responsibility  to  Algorithm  6  which  is  the  modified  version  of  Algorithm  4. 
The  primary  modifications  that  are  present  in  Algorithm  6  are  based  on  the  fact  that  the 
determination  of  the  values  for  H+  and  H~  are  the  key  to  the  problem  as  they  determine 
which  hexads  will  be  created.  If  an  algorithm  can  be  devised  that  can  determine  the 
identifying  history  H  of  the  hexad  that  represents  (as  an  equivalent)  some  other  history  i 
then  the  optimal  set  can  be  produced  directly. 

Algorithm  7  performs  this  function.  It  is  based  on  the  observation  that,  when  equiv¬ 
alent  pairs  of  hexads  (hi,/i2)  occur  in  Sa^  (where  the  convention  of  listing  pairs  in 
descending  order  of  identifying  histories  is  followed),  the  member  of  the  pair  hi  with  the 
largest  identifying  history  h\.H  has  a  component  hi.H+  whose  value  is  not  ★  and  is  smaller 
than  the  identifying  history  (i.e.  hi.i7+  ^  *  and  hi.if+  <  h\.H)  and  that  the  identifying 
history  h2.H  of  the  other  member  of  the  pair  is  equal  to  the  integer  portion  of  h\.H+/2. 
Therefore,  in  order  to  determine  the  smallest  equivalent  history  to  some  history  i,  Algo¬ 
rithm  7  first  determines  whether  i  shifts  to  the  accepting  state  ★  after  a  positive  result  in 
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Algorithm  6  Add  a  hexad  with  H  =  i  to  the  optimal  set  S'M  N 

create  new  hexad  h 

h.H<=i 

h.P  <=  0 

if  #  (z)  —  (M  —  1)  then 
kJB+  4=  * 
else 

n+  <=  ((z  x  2)  mod  2N~l  4- 1) 

h-H+  <=  (Smallest  equivalent  history  to  n+  in  an  M  of  N  automaton)  -  Algorithm  7 
end  if 

n~  ((z  x  2)  mod  2Ar_1) 

h.H  (Smallest  equivalent  history  to  n~  in  an  M  of  N  automaton)  —  Algorithm  7 
$M,JV  S'm}N  ^  {^} 

if  a  hexad  with  H  —  h.H+  does  not  exist  in  S;M  N  then 

Add  a  hexad  with  H  =  h,H+  to  the  optimal  set  S'M  N  -  Algorithm  6 
end  if 

if  a  hexad  with  H  =  h.H~  does  not  exist  in  S'M  N  then 
Add  a  hexad  with  H  =  h.H~  to  the  optimal  set  SrM  N  -  Algorithm  6 
end  if 


Algorithm  7  Smallest  equivalent  history  to  i  in  an  M  of  N  automaton 
if  #  (0  =  (M  —  1)  then 
return  i 
else 

n  <=  ((z  x  2)  mod  27V~1  +  l) 

n'  <=  (Smallest  equivalent  history  to  n  in  an  M  of  N  automaton)  -  Algorithm  7 
if  n1  <  %  then 

return  the  integer  portion  of  n1/ 2 
else 

return  z 
end  if 
end  if 


which  case  it  is  the  smallest  equivalent  state  to  itself.  If  this  is  not  the  case  then  the  history 
n  of  the  state  to  which  z  would  shift  on  a  successful  trial  in  the  non- optimal  automaton 
is  determined  and  then  a  recursive  call  to  Algorithm  7  is  made  to  determine  the  smallest 
equivalent  history  nf  to  n  which  is  the  state  which  i  will  shift  to  in  the  optimal  case.  This 
recursive  call  provides  the  equivalent  to  the  exhaustive  iterations  in  the  method  described 
above.  The  value  of  vl  is  compared  to  z  and  if  it  is  smaller  then  i  will  have  an  equivalent 
state  in  the  optimal  set  of  states  and  the  history  of  that  equivalent  state  is  returned. 


6  Example  of  Application 


Consider  the  case  where  a  surface  ship  equipped  with  an  active  sonar  set  (where  the 
sonar  operator  can  be  considered  a  3  out  of  5  sliding  window  detector)  is  approaching 
a  submarine  and  it  is  required  to  determine  the  probability  that  the  ship  detects  the 
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submarine  before  some  given  range  (for  example,  see  ref.  4).  In  order  to  solve  this  problem 
using  the  methodology  described  in  this  paper  the  following  information  is  required: 

•  the  probability  of  detecting  the  submarine  on  a  single  sonar  ping  as  a  function  of 
the  range  between  the  ship  and  the  submarine, 

•  the  time  interval  between  successive  sonar  pings,  and 

•  the  positions  and  velocities  of  the  ship  and  submarine. 

For  the  purposes  of  this  example,  the  probability  of  detection  on  a  single  ping  (as  a  function 
of  range)  will  be  that  illustrated  at  Figure  3.  This  is  a  fictional  and  stylised  version  of  a 
situation  often  encountered.  There  are  two  range  bands  where  the  probability  of  detection 


Figure  3:  The  probability ,  p,  of  detection  on  a  single  ping  as  a  function  of  range 
is  significant: 

•  a  convergence  zone  between  (approximately)  27  and  37  nautical  miles  (nmi),  and 

•  a  direct  propagation  zone  inside  (approximately)  9  nmi. 

Outside  these  two  zones  the  probability  of  detection  is  roughly  equivalent  to  the  false 
alarm  probability  which  has  been  set  to  (an  unrealistically  high)  0.05  for  the  purposes 
of  illustration.  In  order  to  maintain  a  (convenient)  linear  relationship  between  time  and 
range  it  will  be  further  assumed  that  the  ship  and  submarine  are  on  a  collision  course  with 
a  closing  rate  of  15  knots.  In  order  to  take  advantage  of  the  possibility  of  convergence 
range  detections,  the  time  interval  between  pings  will  be  set  to  that  which  provides  an 
unambiguous  maximum  range  of  40  nmi  (approximately  99.4  seconds)  with  the  first  ping 
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Figure  4-  The  probabilities ,  p*?  that  satisfaction  of  the  3  blips  in  5  pings  criterion  is  met 
at  a  given  range  and ?  P*7  that  satisfaction  of  the  criterion  occurs  by  a  given  range 

occurring  at  40  nmi.  The  probability  (P*)  that  the  3  blips  out  of  5  pings  criterion  is  met  by 
a  given  range  (as  a  function  of  that  range)  is  shown  at  Figure  4  along  with  the  probability 
(p*)  that  the  satisfaction  of  the  criterion  occurs  at  a  given  range.  For  the  convenience 
of  the  reader  the  samples  of  the  blip/ping  ratio  used  are  also  shown  on  the  diagram  as 
individual  points.  It  should  be  noted  that  p*  is  bimodal. 


7  Conclusion 

The  algorithms  described  in  this  paper  for  evaluating  M  out  of  N  sliding  window 
detectors  are  simple  and  efficient.  Their  implementation  in  most  computer  languages 
should  be  simple  and  straightforward.  Applications  include  modelling  of  naval  surveillance 
radar  and  sonar  systems. 
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